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Abstract: We consider probe Dp-branes, p = 3,5, 7, in global AdS^ x S^, rotating 
along an internal direction in the S^. These are dual to strongly interacting A/" = 4 
SYM on 5*^ with massless fundamental flavours, in the presence of an R-symmetry 
chemical potential for flavour fermions. For massless, "AdS-flUing" Dp-brane embed- 
dings at zero temperature, we flnd an inflnite set of threshold values of the chemical 
potential at which instabilities are triggered. The onset of instability is always pre- 
ceded by metastability of the zero density state. From the onset values of the chemical 
potential, we infer that unstable directions favour a homogeneous condensate of a 
bilinear made from fermion harmonics, or Cooper pairing. We conflrm this picture 
both numerically and analytically. The linearized analysis showing the appearance 
of instabilities involves a charged scalar in global AdS space coupled to a (large) 
constant background gauge potential. The resulting frequency space correlator of 
the fermion bilinear at strong coupling displays poles in the upper half plane. In 
contrast, the correlator at zero coupling exhibits Pauli blocking due to occupation 
of states below the Fermi level, but no instabilities. The end-point of the strong 
coupling instability is not visible in our setup. 
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1. Introduction and summary 



The behaviour of strongly coupled systems under extreme densities and temperatures 
is important in a variety of physical situations. Gauge/gravity dualities [jl], H 
present us with an approach towards tackling the physics of such situations, both 
for pure gauge theories |^ and gauge theories with quark-like matter H, ||, 
1^, 0. Holographic duals of strongly interacting fermionic systems at finite density 
are of considerable interest in the context of superconducting/superfluid systems 



np. Such dual descriptions can also provide the necessary window into the workings 
of strongly interacting quark matter. This is particularly relevant for QCD where 
phase transitions are conjectured to occur precisely at densities where the coupling 
is strong, rendering perturbative methods ineffective, and where lattice approaches 
need to surmount the technical obstacle posed by the sign problem. 

In this note we explore properties of one such setup involving "probe" Dp-branes 
{p = 3, 5, 7) in global AdS^ x S^, with the branes chosen to spin along a compact 
internal direction in the 5*^. The dual conformal field theories live on spheres. Im- 
parting a rotation to the Dp-branes in the internal directions allows for a chemical 
potential (conjugate to a U{1) R-charge) for massless fermion flavour fields trans- 
forming in the fundamental representation of the gauge group in the boundary field 
theory. 

One of our original motivations for examining these systems in global AdS space- 
time was to make contact with recent studies of the phase diagram of SU{N) gauge 
theories in finite volume with fermion flavours in the fundamental representation 
[ P^ , |T3|. However, aspects of our setup and the resulting physical picture are pos- 
sibly also of relevance to "top-down" approaches towards holographic superconduc- 
tor /sup erfiuid models involving probe branes (e.g. [l^, |T6|, |l^ [18], |9]). 

The works of |l3l found, using a weak coupling analytical approach at large- 



A^, and small volume lattice simulations of QCD, an exotic phase structure of SU{N) 
gauge theories as a function of quark/fermion number chemical potential. It was 
shown that at large- with a large number of flavours, the systems exhibit an infi- 
nite sequence of third-order Gross- Witten transitions, across which, fermion number 
increases in discrete steps. A jump in fermion number occurs when the chemical 
potential crosses the energy of an elementary fermion harmonic on the sphere, and 
the corresponding energy level gets occupied. This jump is accompanied by a spike 
in the Polyakov loop order parameter and a pair of third-order transitions ^. 

It is natural to ask whether the physics of fermions seen above, can be observed 
in a setting that has a strong coupling description within a gravity dual framework. 
The probe brane configurations we consider present us with such a setup. However, 



^Lattice studies of QCD at finite volume showed tlie same qualitative jumps in quark 
number correlated with spikes in the Polyakov loop. However, these studies also found certain 
qualitative differences including the appearance of metastable configurations. 
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while |T^, 13 1 treated theories with Nf and N comparable, we will only explore the 



theory in the 't Hooft large- iV limit with Nf/N — )■ 0. 

The main result of our study is that the strongly coupled system in finite volume 
does not see phenomena associated to the energy levels of individual, elementary 
fermion states. Instead, the system is unstable to the condensation of compos- 
ite scalars made from bilinears of fermion harmonics on the sphere. The transi- 
tions associated to this occur when the chemical potential times the R-charge ap- 
proaches twice the energy of a fermion harmonic, eventually triggering an instability 
towards fermion bilinear condensation. The strong coupling analysis shows that a 
large enough R-symmetry chemical potential acts as a negative mass squared for the 
fermion bilinears. As expected, the instabilities persist in the infinite volume limit, 
and potentially imply a phase with spontaneously broken U{1) R-symmetry in the 
presence of a chemical potential. A complete description of such a symmetry break- 
ing ground state would provide a holographic model of a superfiuid/superconductor. 
However, within our probe prane ansatz, the instabilities appear to be of a runaway 
nature and the end-point solution is not visible in the setup. Investigations of spin- 
ning fiavour branes in the Poincare' patch of AdS^ were first carried out in [^, while 
a different way of inducing R-symmetry breaking in the spinning brane setup, via a 



probe brane magnetic field, was investigated in [14 



The spinning probe Dp-branes considered here, describe strongly coupled Af = 4 
supersymmetric SU{N) Yang-Mills theory on the three-sphere, coupled to flavour 
multiplets in a way that preserves eight supersymmetries {JV = 2 in four dimensions). 
When the flavours are massless, the theories possess a global R-symmetry. The 
angular velocity of the probe D-branes along a compact internal direction acts as a 
chemical potential for this R-symmetry under which fermions in the flavour multiplets 
are charged. In the strict limit N ^ oo and Nf flxed, which is the probe limit for the 
Dp-branes, the R-symmetry is non-anomalous and is a conserved quantum number. 
Furthermore, since the flavours do not "backreact" on the pure gauge dynamics in 
this limit, we may consistently choose all the R-charge to be concentrated in the 
flavour sector alone. The "leaking" of the charge to adjoint matter flelds of A/" = 4 
SYM is suppressed in the large- limit. In the gravity picture this means that we 
can have the probes spinning along the internal directions, whilst the supergravity 
background is not boosted along the same internal directions. 

In addition, in the probe limit, the corresponding fleld theories remain confor- 
mally invariant. The D7-brane probe yields a four dimensional M = 2 SUSY gauge 
theory, whilst the D5- and D3-brane probes give rise to 2+1 and 1+1 dimensional 
"defect" CFTs respectively ^T], ||. In the global AdS description, the defect CFT's 
are localized on an equatorial 5*^ and 5*^ of the boundary three-sphere, and the 
degrees of freedom localized on the defects interact with the A/" = 4 theory on 5*^. 

The strong coupling analysis of the probe Dp-brane action reveals that the mass- 
less, zero density embedding encounters an inflnite sequence of instabilities with 
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Figure 1: Left : The free energy as a function of chemical potential /x for the D7 system. 
The zero density state (solid black) is locally stable as is increased from zero. A new 
unstable finite density saddle point (dashed blue) with chiral condensate appears when fi 
approaches The two merge at fiR = | where an instability develops. Right: The 

situation is schematically depicted in terms of an effective potential (see also Fig. 8). 



increasing chemical potential /i, for the U{1) R-symmetry. In particular the insta- 
bilities favour the condensation of fermion bilinears with R-charge Qr = +2, made 
up of harmonics of elementary fermions on the sphere. For the D7-branes we find 
instabilities at 

fiQR = {2i+l)R-\ £=1,2,3,... (1.1) 
while for the defect CFT's on S"^ and S^, the respective onset values are 

fxQB = 2iR^^ and ^Qr = {2i - 1) R-\ £=1,2,3,.... (1.2) 

The R-charge of the bilinears is Qr = 2, and in each case the critical value of fiQR, 
is twice the energy of an elementary fermion harmonic (on S^, S*^ and S^). 

The solutions to the non-linear equations of motion of the DBI action, which 
we analyze numerically, display a somewhat intricate structure in the vicinity of 
each of these instabilities. The zero density state, which starts off being stable, 
actually becomes metastable before the onset of the first instability, accompanied by 
the appearance of a new thermodynamically unstable solution with finite density and 
fermion bilinear condensate. This latter saddle point eventually merges with the zero 
density solution at the critical value of the chemical potential. Beyond this value the 
zero density solution is simply unstable. This phenomenon is summarized in Fig. 1 
and in Fig. 8. 

Analytical solutions of the linearized equations of motion for the Dp-branes show 
the emergence of non-trivial embeddings (at linear order) only at the critical values 
of the chemical potential. This confirms the result of the numerical analysis which 
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points towards a merger of two saddle points (a zero density trivial solution and a 
non-trivial finite density solution) at the critical values of the chemical potential. The 
linearized equations which lead to non-trivial massless embeddings, are precisely the 
equations of motion of a charged scalar in global AdS spacetime, coupled to a constant 
background U{1) gauge potential. For large enough, special values of this background 
gauge field (the chemical potential), the equations allow for charged fermion bilinear 
condensates or "scalar hair". In this respect the systems are somewhat similar to 
holographic superconductor models . 



The linearized approach also allows the computation of the two-point function 
for the fermion bilinear which displays an infinite set of discrete poles in the complex 
frequency plane. At = 0, the locations of the poles are at twice the energies of the 
discrete fermion harmonics. With a non-zero chemical potential, the poles migrate 
and eventually move to the upper half plane, signalhng instabilities when onset values 
of the chemical potential are crossed. 

While the detailed features associated to the fermion bilinear instability follow 
from our anstaz for the probe D-brane action (which only incorporates the so-called 
"slipping mode" degree of freedom), the end-point of these instabilities does not seem 
to be captured by this action. This suggests that new effects must be included to 
obtain a complete description. 

The paper is organized as follows. For the purpose of completeness, in Section 2, 
we review basic field theory features of the probe D7-brane system, paying attention 
to the introduction of an R-charge chemical potential for flavours alone. In Section 
3, we explain weak coupling features of the theory on a sphere with non- vanishing 
chemical potential. The probe D7-brane analysis is done in Section 4, with numerical 
results for the phase diagram and free energies of solutions. Analytical results in the 
linearized approximation and strong-coupling Green's functions are derived in Section 
5. We discuss the rotating D5- and D3'-branes, related to D3/D5 and D3/D3' systems 
in Section 6, and conclude with speculations on the end-point of the instabilities in 
Section 7. Various details of weak coupling calculations of the phase structure on 
5''^, Green's functions at finite density, and technical features of the D7-brane ansatz 
are presented in an Appendix. 



2. Adding flavours to A/" = 4 SYM on 5^ x R 

It is well known that the simplest way of incorporating flavour fields into the AdS/CFT 
framework is via probe Dp-branes in the bulk type IIB supergravity background dual 
to, for example. A/" = 4 SYM with SU{N) gauge group 0. The "probe" approxi- 
mation is valid when the the number of flavour multiplets Nf is fixed in the 't Hooft 
large- limit. 

The presence of flavour fields spoils conformality of the A/" = 4 theory since the 
gauge coupling now runs and has a putative Landau pole in the ultraviolet (UV). 
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In the limit Nf/N — however, the loop contributions from flavour modes are sup- 
pressed and the beta function for the 't Hooft coupling, A = ^fy^iV, is proportional 
to ^ and therefore vanishes in the strict large N limit. We may continue to view the 
theory as being conformal, provided the hypermultiplets have vanishing bare masses. 
It is natural for CFTs in four dimensions to be viewed as fleld theories formulated 
on a spatial three-sphere, a picture that follows from radial quantization on M^. 
The radius of the three-sphere provides a new scale, allowing for the possibility of 
non-trivial thermodynamics as a function of temperature and/or chemical potentials 

MM- 



There are three distinct cases to be considered, corresponding to probe D7-, 
D5- and D3-branes. In all cases the resulting system preserves 8 supersymmetries. 
The D7-brane probes describe the dynamics of M = 2 SUSY gauge theory with Nf 
hypermultiplets, while the dynamics of the two other kinds of probe branes yields 



lower dimensional defect conformal fleld theories pT] , p2| . We will begin with the 
study of the D7-brane probes. 

2.1 D7-brane probes: Matter content and symmetries 

The matter content and the interactions of the fleld theory are most easily summa- 
rized in the language of = 1 superflelds on M^'^. Each M = 2 hypermultiplet 
consists of two M = 1 chiral superflelds (Q, Q) transforming in the [N, N) repre- 
sentation of the SU{N) gauge group. The superflelds in the hypermultiplet, each 
consist of a squark and a quark Q = (g, ipg) and Q = (g, ipg). The quarks ipg and 
are left- and right-handed Weyl fermions respectively. 

With Nf hypermultiplets {Q\ Qi) and using the M = 1 SUSY notation, the 
superpotential of the theory with flavours on M?'^, can be written as 



9ym 



5^(y2 4$3Q* + m4Q*) +y2Tr($3[$i,$2]) ) . (2.1) 



The six adjoint scalars of the A/" = 4 theory, and their superpartners are packaged 
into the three (complex) chiral multiplets $i,$2 and $3. Upon coupling to the 
M = 2 hypermultiplets as in the above superpotential, $3 naturally gets identifled 
with the scalar component of the M = 2 vector multiplet, while $12 together make 
up an A/" = 2 adjoint hypermultiplet. 

Although one may allow for bare masses for the hypermultiplets, we will only 
be interested in the theory with massless quarks. In this case the theory has the 
following global symmetry group 

Gf ^ U{1)b X SU{Nf) X SU{2)r x U{1)r x SU{2)^. (2.2) 

The multiplets transform in the fundamental representation of the flavor symme- 
try U{Nf) ~ U{1)b X SU{Nf), while Qi transform in the conjugate representation. 
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The SU{2)r X U{1)r factor is the R-symmetry group of the Af = 2 gauge theory. 
The SU{2)ji subgroup acts on the doublets and leaves all other fields in- 

variant. We will be interested in the U{1) subgroup of the R-symmetry. Under a 
U{l)ji transformation, the scalar component of $3 has charge +2, whilst the squarks 
in the hypermultiplets are uncharged. Importantly, the fermions ijj^^ and ^^^^ each 
carry R-charges -|-1 and —1 respectively. The ?7(1)_r rotation is therefore an axial 
symmetry acting on the fermions ^. 

The U{l)ii symmetry is anomalous for a fixed number of colours and flavours. 
However, M = 2 supersymmetry relates this anomaly to the beta function for the 
gauge theory. Like the beta function, the U{1)r anomaly is also parametrically 
suppressed in the large N limit with Nf fixed. It is thus consistent to treat the 
U{1)r as a symmetry of the theory in the 'tHooft large- limit. 

2.2 The theory on 5"^ x M 

When the N = 2 SUSY gauge theory described above is formulated on a spatial 
three-sphere, conformal invariance requires that we add to the Lagrangian of the 
theory, a coupling of the scalar fields to the Ricci scalar of S^: 



3 Nf 
a=l 1=1 



:2.3) 



Here R is the radius of the boundary three-sphere. The conserved current associated 
to the U{1)r symmetry is a sum of two parts: a contribution J a from the charged 
adjoint matter of the Af = 4 sector and another piece Jp from the R-charged funda- 
mental fermions. 

= J'a + J'f (2.4) 



= i Tr (2 <I>;^^''$3 + ^^a"^ + A^cr'^A 



Chemical potentials for global symmetries can be turned on in the Lagrangian for- 
mulation by first imagining that the corresponding global symmetry is gauged and 
then introducing a background value for the time component of that gauge field. For 
the U{l)ji symmetry above, a chemical potential will shift the Lagrangian as, 

+ j°)-/i2|$3|2. (2.5) 



We take the fermions i/'g and -01^ ,j to be left and right handed chiral fermions respectively. This 



is consistent with the Yukawa coupling implied by the superpotential (2.1) 
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As is well known |^6|, ^ for yU > and in the absence of fundamental flavours, 
the chemical potential induces a runaway instability in the potential energy for the 
adjoint scalar $3. However, this issue will not be a problem for us as we explain 
below. 

As we will work in the probe limit for flavour branes, the number of flavour 
hypermultiplets is small, Nf <^ N, so that the flavours do not "backreact" on the 
jV = 4 gauge theory. We then introduce a chemical potential for the f/(l)_R charge 
in the flavour sector alone. Strictly speaking, this will not result in an equilibrium 
configuration since the R-charge carried by the flavour degrees of freedom will even- 
tually leak out to the R-charged adjoint modes of the full theory ^. However, in the 
't Hooft large limit with Nf fixed, this effect is suppressed. 



Figure 2: Typical graphs that contribute to a one-point function for Jj^, the R-charge 
density in the adjoint sector, when a chemical potential for the flavour R-charge density 
is introduced. ^3 and ip are the complex adjoint scalar and its superpartner, which are 
both charged under U{1)r, while {q,ipq) represent flavour squarks and quarks respectively. 

As shown in Fig. @, a chemical potential conjugate to the flavour R-charge 
density Jp also implies a non-zero one-point function for R-charge density for adjoint 
matter. These effects are suppressed by flavour loops in the 'tHooft large- A^ limit. 
Therefore in this limit, it will be consistent for us to introduce an R-charge chemical 
potential in the flavour sector alone, 

C^C-fiJ^. (2.6) 



3. Weak coupling interlude 

For the weakly coupled theory on 5''^, the energies of all the perturbative modes are 
quantized and there are no zero modes (at zero temperature). The putative zero 



^The spinning D-brane setup has been applied to understand aspects of this non-equiUbrium 



system with massive hypermultiplets in flat space [29 . In this case the rotation of the branes gives 
rise to a time-dependent phase for the mass. 
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modes for scalars are lifted by the conformal coupling of the scalars to the curvature 
of the sphere. The eigenf unctions of the spinor Laplacion on have the eigenvalues 
and degeneracies (see e.g. [ pO|] ) 

e^=[i+l)R-\ de = 2£{£+l), £=1,2,... (3.1) 

At large- A^, for a non-anomalous R-symmetry, the chemical potential for the U{1)r 
symmetry in the flavour sector is conjugate to the chiral fermion number, which 
is the number of left-handed Weyl fermions minus the number of right-handed 
fermions Nr. At zero temperature, as fi is increased past each of the energy levels 
( p.l|) , the quantity Nl — Nr will undergo a discrete jump, 

A{NL-NR) = deNfN, (3.2) 

as the energy levels below ^ get filled. 




Figure 3: Left : The low T, phase structure at zero coupling and large with fixed Nj/N 
on S^. Right: Behaviour of — Nji and the Polyakov loop with increasing fi at low T. 



3.1 Weak (zero) coupling with Nj/N fixed 

The effect is particularly intricate in the Veneziano limit (assuming for the present 
purposes that the axial U{1) symmetry is anomaly free) on the three-sphere. At non- 
zero temperatures, the Euclidean theory on x 5*^ reduces to a unitary matrix model 
for the Polyakov loop matrix U = expi Aq. This effective theory is obtained by 
integrating out all massive Kaluza-Klein modes on (e.g |]25|). As is well known 
[p5[] from this analysis, the free theory without fundamental matter experiences a 
Hagedorn/deconfinement transition at large- A^. In particular, at low temperatures 
the theory is in a confined phase with j^{TrU) = 0, while beyond the deconfinement 
transition this order parameter acquires a non-zero expectation value. 

Incorporating a large number Nf of fundamental matter fields with Nj/N fixed, 
changes the low temperature picture, especially when a chemical potential for the 



fermions is introduced. The situations with baryon number chemical potential [12 
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and f/(l)i? chemical potential (see Appendix ^ for the fermions are quite similar, 
although the respective technical details differ. 

From the analysis in Appendix 0, we see that as the flavour R-charge chemical 
potential /i is increased towards the energy level e£ of a fermion mode, the theory 



experiences a third-order Gross- Witten transition pl| at 



~ ee-T\n(^^d,j . (3.3) 

Here, the low temperature approximation implies that eg ^ T. At this transition, 
the large- distribution of the eigenvalues of the unitary matrix U changes from 
an ungapped to a gapped phase. The distribution, however, quickly reverts to its 
ungapped phase as is increased shghtly past eg, at 

fil ~ £, + Tln(f rf,) (3.4) 

which represents yet another Gross- Witten transition. This pair of transitions is 
accompanied by two qualitatively interesting effects. The first is a change in the 
Polyakov loop order parameter -^{TiU) which is small in the ungapped phases, but 
becomes large in the gapped phase. The second physical effect is the chiral fermion 
number which jumps across the two transitions so that A{Nl — A"^) = d^NNf. 
These phenomena are depicted in Figure (^. The work of has shown that these 
qualitative effects (with a baryon number chemical potential) persist in QCD (A^ = 3) 
in a small volume. 

3.2 Green's function for ipqipq at weak coupling 

The study of the large- A^ theory at strong coupling and with Nj/N fixed is beyond 
the scope of this paper. The main point of the weak coupling picture described 
above is to indicate the various effects associated to threshold values of the chemical 
potential. As we will see below, the elementary fermionic states are not visible 
at strong coupling; instead, the fluctuations of probe D-branes at strong coupling 
comprise of composite gauge invariant operators made up of fundamental fermions 
and their superpartners. The simplest amongst these is the fermion bilinear 

Nf 

(3.5) 

i=l 

This is a dimension three operator at weak coupling. In the context of the D3- 
D7 system its conformal dimension is protected by its R-charge 0. However, the 
normalization of its correlator can and does receive quantum corrections as a function 
of A = Qym^ ■ At strong coupling, in the presence of an R-charge chemical potential, 
we will be interested in computing the two-point function of this operator on . It 
will be extremely interesting to contrast this with perturbative results. The frequency 
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Figure 4: Two-point correlator of ipqipq in the free theory, in frequency space. 



space retarded correlator of O, at zero external spatial momentum is, 

Qnico) = ^J^°^dte^-'J ^ , t) , 0(0)]) . (3.6) 

Using the finite density propagators ( p.5| ) for the internal fermions in Fig.(|D, we 
obtain a general form for this correlator at zero temperature and non-zero chemical 
potential /x > 0, 



For the retarded correlator we need to choose an "ze" prescription where the poles 
in the correlator are taken to be slightly below the real axis. The simple poles in 
u) represent on-shell intermediate modes with vanishing external spatial momentum. 
The step function B(/i — e^) is the Fermi-Dirac distribution evaluated at zero temper- 
ature, and exhibits the "Pauli-blocking" of occupied levels below the energy /j,. Using 
the energy levels ei and degeneracies for S"^, the Green's function can be written in 
terms of digamma functions (up to constants and "contact" terms proportional to 



n 



-2^ (I + n + looR) - 2^- (I + n - ^ioR)] 



Here n is the largest integer less than or equal to fiR — |. In the perturbative 
limit, the chemical potential does not introduce new poles in the correlator Q{u). In 
particular, as expected there are no instabilities of any kind. In the limit of large 
R, we should recover the correlator in flat space, after incorporating the correct ie 
prescription. 



167r^ 



\xi{u?R^) + In ( ^^^^j^) + 2m sgn(u;) Q{y? - 4/^^) 



(3.9) 
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using il>[x) — )■ Inx for large x. The correlator has branch points at w = ±2yU and 
an imaginary part for u"^ > Afi"^. This is because states with energy below \^\ 
are occupied and therefore Pauli-blocked. Only when the external energy \u\ is 
larger than 2|/i|, can the intermediate states be on-shell and yield an imaginary part 
for the Green's function. When /i = 0, the correlator has the expected behaviour 
~ u'^lnu) for a dimension three operator. In particular, at /i = we expect the 
two-point function (including when R is finite) to have the same form for all values 
of the 't Hooft coupling , although the normalization can and does receive quantum 
corrections. We now turn to the strong coupling description of the system in the 
probe limit Nf <^ A^. 



4. Probe D7-branes in Global AdS 

In the large- limit and at strong 't Hooft coupling A = gyM^ ^ ^i flavour 
modes correspond to low-energy excitations of probe D7-branes in the AdS^ x ge- 
ometry [^. The flavour D7-branes wrap an C while filling the AdS^ directions. 
Such embeddings of the D7-brane are characterized by the "slipping angle" 9, which 
determines the size of the three-sphere wrapped by the D7-brane and is defined by 
the following paramerization of the metric 

dnl = de^ + siTi^ed(t)^ + 008^6 dnl, o<^<-. (4.1) 

A static probe D7-brane wrapping this breaks the 5*0(6) R-symmetry of A/" = 4 
theory to an SU{2)r x U{1)r subgroup. As discussed above, the U{1)r is non- 
anomalous in the theory with ^ — )■ and is the isometry generated by shifts of 
the angle 0. We will consider probe branes spinning in the U{1)r isometry direction 
parametrized by the angle (p. The associated angular velocity is conjugate to the 
U{1)r charge carried by such states. 

When embedded in the global AdS spacetime, the D7-brane world-volume is 
topologically S*^ x S"^ x M x R^. The first 5*^ sits inside the global AdS^ geometry 
whose conformal boundary is an . The D7-brane wraps a second S*^ C 5*^, specified 
by the polar angle 9 in Eq. (^4.1|) . We write the global AdS^ metric in Fefferman- 
Graham coordinates 



ds' = -^i^All dt' + dnl + %, (4.2) 



with < z < 2. The conformal boundary is approached as z — )■ 0. We have set the 
AdS radius to unity. In these units and the coordinate system employed above, the 
radius of the boundary three-sphere is R = 1. 
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4.1 Embedding Ansatz 

We can take the D7-brane to spin with constant angular velocity in the direction, 
while at the same time allowing for a radial dependence of its angular coordinate on 
the 0, ^ 

(f)= fit+ g{z). (4.3) 

It turns out (Appendix C) that for massless flavours the radial dependence must be 
taken to vanish. So we will only consider the configurations with g{z) = 0, as we are 
mainly interested in the m = theory when the U{1)r can be viewed as a symmetry 
(at least in the large- A limit). The angular velocity jl is related to the R-charge 
chemical potential in field theory, as 

fL = 2fi. (4.4) 

The factor of two appears because the two fields excited on the probe brane, namely 
6 and are related to the magnitude and phase, respectively, of the fermion bilinear 
made from flavour fields. This object has R-charge +2. An R-charge chemical 
potential for a fundamental fermion flavour can also be viewed as a time dependent 
phase exp{ifit) for the fermion fields. Correspondingly, the bilinears have a time 
dependent phase exp(2i/it). 

Setting g{z) = for massless flavours, the induced metric on the D7-brane is 



D7 



(4.5) 



_ I 1^+4^ ) _ 4^2 gi^2 ^ ] ^ (z^\zf + 1) ^ + ^ i ' + COS^ Odnl 

\ Z^ J z^ 

The D7-brane probe dynamics is then controlled by the DBI action 

^D7= 271^07 dtdz cos^e{z)- ^ {1 + z^O' {z)^){{l + ^ z^)^ - A/j.^ z^ sin^ 6 

(4.6) 

where Ad7 depends on the D-brane tension T^-j and the number of flavours Nf, as 

Afnj ^ 2n' Nf T^, = ^^,NfN . (4.7) 



The additional normalization factor of 27r^ in Eq. ( |4.6| ), can be traced to the volume 
of the boundary three-sphere, whose radius is 1 in our units. 

Classical solutions to the equations of motion resulting from varying Eq. (|4.6| ) 
have the following general behaviour near the AdS boundary, 

Oiz) 1,^0 = ^(0) z + 6(^2) + \ ^(0) (1 - 4/i2) zHYiz + .... (4.8) 

From the standard AdS/CFT dictionary this is the behaviour expected for a 
supergravity mode dual to a dimension 3 operator in the field theory. It is interesting 
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to note that in global AdS, there is an extra logarithmic term in the asymptotic 
expansion which explicitly depends on the chemical potential. This will give rise 
to a chemical potential dependent counterterm for the probe brane action below. 
The operator must be made from flavour fields and must transform under the global 
U{1)r symmetry. It is natural then, to associate the slipping angle 6 with the quark 
bilinear^ 

= J2^,^^l■ (4-9) 

i=l 

The leading coefficient ^(o) , when non-zero, results in a bare quark (hypermultiplet) 
mass which explicitly breaks the t/(l)_R symmetry. For vanishing hypermultiplet 
mass, the coefficient of the term directly yields the quark condensate (O). A non- 
zero quark condensate breaks the f/(l)_R symmetry spontaneously (assuming that the 
quark mass is zero). The precise dictionary between the gauge theory parameters 
and supergravity modes is quoted below. 

4.2 Divergences and renormalization 

The D7-brane probe action is formally divergent and needs to be regulated and 
renormalized. The regularization of these divergences is well understood within the 
general framework of the holographic-RG approach applied to probe D-branes [p4|. 



35| , |36| . It is instructive to first look at all the divergent counterterms implied by 



the near-boundary asymptotics (4.13) in the D7-brane action, and then compare 



with the results of the holographic RG method. By cutting off the z integration at 
2; = e ^ 1, we find the following divergent terms 



S^, ^ 2n'Afj,r I dt ^-1-^ + _L + !| + ^2^^ (1 _ 4/,2) Ine + . . . j , (4.10) 

Note that the chemical potential /j, contributes to the logarithmically divergent coun- 
terterm. However, when ^(o) = and the hypermultiplets are massless, the diver- 
gences are /^-independent. This is consistent with the general expectation that much 
like a finite temperature, a chemical potential should not introduce new divergences 
in the theory. For non-zero hypermultiplet masses, U{1)r is not a good symme- 
try and the angular velocity /x leads to a time-dependent mass and also makes an 
appearance in the UV divergent terms. 

The divergences are removed by the addition of local counterterms at z = e to 
yield the regularized action 



r ^ 

S'^^ = S^j + MYi7 / rfl^3 dt Lk (4.11 

k=i 



■*The precise mapping relates 9 to the operator i^qii^^ + <z|$3 9* + <7i$3 <7^' 
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Denoting the boundary metric of the cut off AdS space as ^fj^ydx^dx" = -^{—dt^ + 
c/fig), we can write these counterterms in terms of the Ricci scalar TZ of the boundary 

at z 



BR ^ 



e as m 
1 



-7 



Lo = 



Lfi = — 



12 



2^ 



-eief (4.12) 



We have used the fact that ^(e) is given by Eq.( 4.13 ) as e — )■ and that TZ = 6e^. 
Notice that the counterterm L5 depends on the chemical potential fi. The form of 
this counterterm can be deduced by first noting that that 6 and can be packaged 
into a complex scalar $ and the counterterms for this scalar can be expressed in the 
form of Eq.(4.2) in [||] ^. 



Having defined the regulated action, the dictionary between the gauge theory 
parameters, namely the quark mass and condensate, and corresponding supergravity 
modes can be worked out easily (having set the AdS radius to unity). Specifically, 
for massless hypermultiplets we obtain 

" (^A) =A/'d7^(2) (4.13) 



m 



0. 



4.3 Solutions 



The DBI equations of motion have a large family of spinning solutions. The solutions 
fall into two topologically distinct varieties. In one class of solutions, the 5*^ C 5*^ 
which is wrapped by the D7-brane probe, shrinks in the interior, before the brane 
reaches the centre of the global AdS space. This type of solution always yields a 
non-zero hypermultiplet mass which can be read off from the asymptotic behaviour 
of 9{z). The second category of solutions are those where the second 5*^ embedded 
in AdS^ shrinks, and the probe brane reaches the origin of AdS. A very interesting 
topology-changing, but continuous phase transition between these two classes of so- 
lutions was discovered in ||36|, |3^. This phase transition was originally found and 
explained at zero temperature, at a critical value of the hypermultiplet mass. It 
can be thought of as a finite volume, meson binding/unbinding transition at strong 
coupling and zero temperature. The second category of solutions - the AdS-filling 
embeddings, include configurations corresponding to massless hypermultiplets. In 
this phase, a maximally separated static quark-antiquark pair is stable against light 
quark pair production. We will be interested primarily in the AdS-filling massless 
embeddings which do not explicitly break the f/(l)i? symmetry. 

4.3.1 Constant Solution 

The simplest solution to the DBI equations of motion following from Eq. ( [4.6|) is 
6{z) = 0. This has m = (O) = 0. The free energy or the action per unit time for the 

would like to thank Andy O' Bannon for pointing this out to me. 
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constant solution is 



D7 



16 



ArD7. (4.14) 



The induced metric on the D7-brane is simply (global) AdS^ x . The charge density 
conjugate to /i 

/■^ dC 

d = M,y, dz— (4.15) 

Jo on 

evaluates to zero for the 6 = Q solution. 

The trivial solution with m = and zero R-charge density {0{z) = 0) exists for 
all values of /i, but as we will see below, it is not always thermodynamically stable. 

4.3.2 Non-constant solutions 

We will now examine in detail the possibility that there may be additional, non- 
constant solutions to the DBI equations of motion with m = ^(o) = 0. If such 
embeddings were to exist, they would correspond to configurations with a finite R- 
charge density and which also spontaneously break the U{1) R-symmetry by the 
formation of a fermion bilinear condensate 6(^2) 7^ 0. 

To find out whether such configurations exist, given the non-linearity of the DBI 
action and equations of motion, we first employ a brute force numerical approach. 
For non-constant solutions which fill AdS space ^, smoothness at the origin of AdS^ 
requires 

e{z) I = T = constant ; 9\z) | = ^ • (4- 16) 

Furthermore, reality of the DBI action ( |4.6| ) imposes the restriction 

sin^(z)|^^^<i- (4.17) 

so that the action is real. At this point we stress that we only focus attention on 
the solutions having = 2/it + g{z) with g{z) = 0. This is because we want to 
isolate the m = embeddings and for these g{z) must necessarily vanish, as argued 
in Appendix C. For large enough /i, we may expect 6{z) to be small everywhere, and 
apply a linearized approximation to the equation of motion. We will return to this 
shortly. 

The full space of non-trivial spinning D7-brane embeddings for general /i, in- 
cludes both the AdS-filling solutions and the solutions which cap off before getting 
to the orgin of AdS^. The latter class of embeddings are characterized by the fact 
that the C wrapped by them shrinks smoothly before the brane gets to the 
centre of AdS^. 

TT 

6{zo) = - , 0'{zo) -> 00 , for some zq < 2. (4.18) 



^Solutions which cap off before getting to the centre of AdS spacetime, always correspond to 
massive flavours, and fall in the category of so-called "Minkowski embeddings" . 
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By numerically solving the DBI equations of motion for different values of /i, we 
extract the mass parameter m and condensate (ipq'ipq) from the asymptotics (|4.13| ). 
In this way we can obtain all possible types of D7-brane embeddings which have the 
slipping angle 9{z) as the only non-trivial mode. The possible space of solutions is 
summarized in Fig.(^ 




Figure 5: The space of known spinning D7-brane embeddings in global AdS^ x S^. The 
known solutions all lie to one side of the curves shown. The figure also shows that near 
fiR ~ |, there is a solution with m = 0. (Note that we have reinstated the radius R of the 
boundary sphere, which was set to 1 by our metric conventions). 



The figure exhibits three "phase boundaries" : 

• The line in blue separates "AdS-filling" solutions from the "capped-off" em- 
beddings. This boundary line is obtained by looking at AdS-filling embeddings 
with 9(z)\_^^^ — i- 7r/2 and 9'{z)\^^^ = 0. These are on the verge of getting 
capped off and form a phase boundary for a continuous topology-changing 
transition of the type found originally in ^ |. 



The red line is the maximum allowed value of the chemical potential fi, below 
which the DBI action still remains real, for the massive capped-off solutions or 
so-called "Minkowski embeddings". This is obtained by noting that a D7-brane 
ending at z = Zq < 2 with 9{zq) = 7r/2 must have (1 + ^^^o) — ^At-Zo for every 
allowed value of mass (from ([4.6|)). 



The third boundary line in black is of interest to us and it represents solu- 
tions that fill AdS^, with the maximum allowed value of /i dictated by re- 
ality of the DBI action ( [4.6|) . These saturate the inequality ( [4.17 ), so that 



2yU = 1/ sm9{z)\z=2- A key feature of this family of solutions is that it inter- 
sects the m = axis at several points close to (but not precisely at) odd integer 
values of fiR, (Fig. 
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Figure 6: Left : A blown-up plot of the masses and chemical potentials for three families of 
AdS-filling solutions. The black line is the boundary line depicted in Fig. 4. New massless 
embeddings appear in the vicinity of /xi? = |, |, | . . .. Right: Closer examination reveals 
massless (m = 0) finite density solutions for a finite range of 2fiR, close to the odd integers. 

Blowing up the region for small m, we see clearly in Fig. (U), that there exist non- 
trivial solutions with m = 0. By plotting the masses for four families of solutions, we 
see that the non-constant, vanishing mass embeddings occur near 2/ii? = 3, 5, 7, . . .. 
It is also clear from the plots that there is a finite range of values of 2/ii? in the 
vicinity of the odd-integer points, for which non-trivial massless solutions arise. For 
example, it seems that such finite density massless solutions to the DBI equations of 
motion appear for 2.9 < 2fiR < 3. 

The four families of solutions, whose masses are displayed in Fig.(|]), were ob- 
tained by integrating the equations of motion outwards from the origin {z = 2) with 
boundary conditions, assuming that 2/i > 1: 

^'(^) L=2 = , eiz = 2)= sin~i n = 1 , 2 , 15 , 35 . 

The values of n were chosen to illustrate the full range of fi in Fig.(^, in the vicinity 
of 2/ii? = 3, for which new m = embeddings exist. The limiting case for which 
the DBI action still remains real, corresponds to n = 1 (solid black curve). The 
families labelled by n = 2, 15 and 35, all satisfy Eq.( [4.17|) and are represented by the 
dotted-black, solid-blue and dotted-red curves, respectively in Fig.(||). 

All the configurations above (including those with m = 0) will necessarily have 
a non- vanishing chiral condensate (ipgipq), indicating that the solutions with m = 0, 
dynamically break the chiral U{l)ji symmetry (see Fig.(|^)). What remains to be 
established is whether these configurations are stable and whether they are thermo- 
dynamically favoured over the constant 6 = solution. 

4.3.3 Free energy and saddle points 

We now turn to the key issue, namely, how the free energies of the non-trivial D7- 
brane embeddings discussed above compare with that of the constant ^ = embed- 
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Figure 7: Left : The normalizable mode 0(2) which is proportional to the chiral condensate 
(ipqipq) when m = 0. These are the values extracted from the n = 1 family of solutions. 
Right: The non-constant massless solution for 2/ii? w 2.95. 



ding. We first note that, as /ii? is dialed from zero, = remains the only solution of 
the DBI equations of motion. A non-constant embedding first appears at ?a 2.9 
(from the n = 1 family above). 

Using the regulated form of the action ( |4.11| ), we numerically evaluate on each 



of the m = solutions, the difference of the regulated actions, 

AS = (27rWD7)-' [S'^'^m = 0,/ii?) - S''%J ■ (4.19) 

In all cases we find that the non-constant embeddings in the window 2.9 < 2fiR < 3 
have a higher action than the 6 = solution. Thus these must either be unstable or 
metastable. In addition, in the limit that 2fiR — )■ 3, the non-constant embeddings 
approach ^ = 0. 

The complete situation is best described by the following procedure leading to 
Fig.(|D. For every non-trivial solution 0{z), obtained with boundary conditions of 
the form ([4.19|) , we define a one-parameter family of configurations which do not 
necessarily solve the DBI equations of motion, 

9{z;s) ^ se{z)[^^^^^^^ 0<s<s,, (4.20) 



s*0{z)\^^^= sin 1 Q 



2/i 

This is a simple rescaling of a given non-trivial solution, where the rescaling parame- 
ter is bounded by the requirement that the DBI action evaluated on the configuration 
be real. It yields a one dimensional slice of the configuration space of 6{z) (with the 
condition that the hypermultiplet mass evaluates to zero at the boundary), along 
which, only the two points s = and s = 1 represent saddle point configurations 
or solutions of the DBI equations of motion for a fixed /i. It is important that s is 
bounded by some value beyond which the DBI action evaluated on these rescaled 
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configurations ceases to be real. Hence the configuration space of 6{z) has a boundary 
for some fixed /x. 

In Fig. (1) we plot the action evaluated on these configurations and find that 
the non-constant solutions are always unstable. In contrast, the zero density, ^ = 
embedding is stable for low /i, but becomes metastable as fi is increased, eventually 
becoming unstable at 2fiR ^ 3. The end-point of the metastability/instability ap- 
pears to lie outside the region of validity of the ansatz for our probe D-brane setup. 
The main result we take away from this analysis is that the ^ = solution becomes 



AS AS 





Figure 8: AF, the difference in free energy (action per unit time) between non-constant 
massless configurations and the zero density state, evaluated along a one-parameter family 
of configurations labelled by s. The point s = represents 6{z) = 0, while s = 1 corresponds 
to the new non-constant solution which first appears when 2fiR > 2.9. The four plots 
illustrate what happens when is increased. At 2fiR = 3, the critical point, the two 
saddles merge and the system has a runaway instability. 

metastable and eventually encounters a runaway instability at 2/ii? = 3. 

The entire picture described in this section is repeated for every odd-integer 
value of 2/ii?. The unstable direction at each of these critical points is towards a 
non-constant massless embedding 6{z). The critical values 2/ii? = 2i + 1, indicate 
that the modes condensing are fermion hilinears which have R-charge 2 and are made 
from elementary fermion harmonics carrying energy (£-|- The mass to charge 

ratio of the bilinears is -|- \)R~^ and governs the critical value of /ii? for the onset 



- 20 - 



of a given instability. Schematically, the operator ipqiipg can be expressed as a sum 
over spherical harmonics, whose homogeneous mode on has the form, 

^ f d'n = (4.21) 

^ e,n 

with ipe^n an eigenmode of the Dirac operator on with mass (£+|)i?~^, £=1,2,... 
and n denotes the azimuthal quantum numbers associated to a given harmonic. 
Therefore, the instabilities are triggered when the chemical potential times the R- 
charge exceeds the mass of the composite charged scalar Oi ~ ^„ V'£,n'^^ 

The unstable directions seen above are also qualitatively distinct from the run- 
away behaviour encountered in A/" = 4 SYM for large enough R-charge pB, pi, E8 



where the effect is due to the presence of flat directions for R-charged scalars. In 
that case, the scalar "Coulomb branch" instability is also visible at weak coupling. 
In the present setup, the (axial) R-charge is carried by the flavour fermions only. In 
the weakly interacting theory, when the chemical potential hits half-integer values, 
the chiral fermion number jumps and there is no instability (see Section(^). In 
contrast, the elementary fermions are not visible in the strongly interacting limit; 
instead condensates of fermion bilinears are preferred when the chemical potential 
approaches the same half-integer points. At strong coupling these values are natu- 
rally interpreted as the mass to charge ratios of fermion bilinear modes. Surprisingly, 
the DBI action, at least within our ansatz, does not appear to include the end-point 
of this instability. Below we speculate on possible resolutions of this. 



5. Analytical results: instabilities 

We now turn to an analytical approach to reveal the origin of the instabilities we 
have encountered above, and in the process gain some understanding of the nature 
of the numerical solutions. The main features of the numerical approach, such as the 
appearance of new albeit unstable classical solutions, are a consequence of the non- 
linearities of the DBI action. However, we expect some quantitative aspects to be 
accessible in a linearized regime, especially for large /i. This can be understood to be 
a consequence of Eq. ( [1.1 7| ) , which requires that for large enough fiR, the amplitude 
of 6{z) is essentially bounded by l/{2fiR). The metastabilities and instabilities we 
saw above appear at an infinite set of values 2fiR = 3, 5, 7, . . ., and so we expect 
some of this physics to be accessible in the linearized regime. 

5.1 Non-constant massless embeddings 

The first issue we would like to understand, is the appearance of a new saddle point 
or solution to the equations of motion with a finite density. When yU ^ 1, sm6{z) < 
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(2/i) ^ <^ 1, and it is then consistent to linearize the DBI equations of motion: 

As we will see more explicitly in the next section when we compute the two-point 
function for fluctuations around the 6' = solution, this is the equation of motion for a 
charged scalar coupled to a constant background gauge field in global AdS spacetime. 
It turns out that the equation can be exactly solved in terms of hypergeometric 
functions. This is somewhat easier to see after the change of variables 

z = 2^, ^(t)^^-^, /(^)^(^_ 1)^/^3^, (5.2) 
following which we obtain a Schrodinger equation for ip{y), 

+ 173 TT,. 73 73^ = 0- (5-3) 



The equation is solved by regularized hypergeometric functions 



Oiz) = C, v^(|i±^ 2F1 (^/i - i, /X + i; 2/x + 1; (i±|j J (5.4) 



where the coefficients are determined by requiring regularity at the origin of AdS 
spacetime at z = 2 or y = 1. Specifically, we are interested in embeddings that fill 
AdS 5 and get to the origin of the space, satisfying 

e{z)\^^^ = T, ^'(^)L=2 = 0- (5-5) 

Imposing these boundary conditions we find 

2fi- 1 7rtan(jr,^) • 2,^+lr(i+^)" 

It is readily checked that the above linear combination yields a real function^ which 
is regular at z = 2. Near the AdS boundary {z — )■ 0), we can then read off the 
coefficients in the asymptotic expansion ( [4.131) . After reinstating the radius R of the 
boundary three sphere which was effectively set to 1 in our conventions, we obtain 



''To check the z — > 2 asymptotics and reahty of the solution, one must choose the correct branch 
of the hypergeometic function which has a branch point at z = 2 (when its argument diverges). 
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the following analytic expressions for the mass m and the condensate (ipqipq) 



8 cos TTfiR 

n 4/i2i?2 _ 1 



^(0) = -T r T-TB^ > (5-7) 



COSTT/ii? 

X [3 - (Vi?2 _ 1) 1^ ^1 _ + ^ (1 + + + - l}] . 

Here is the digamma function which has simple poles at x = —n, for all n G Z. 
Notice that in the expression for 6^(2) above, these poles are precisely cancelled by the 
zeroes of the cosine and 6(^2) is finite for all values of /i. 

Importantly, the set of linearized solutions, parametrized by T, yields massless 
hypermultiplets precisely when 

^0 = ^ 2fiR = 2£+l, £=1,2,3,... (5.8) 

Therefore, non-constant massless embeddings first appear for half-integer values of 
the R-charge chemical potential. This is in agreement with the numerical results 
from the complete DBI action, which exhibit a more intricate structure including 
metastable and unstable regions in the vicinity of the same critical values. It is also 
easily checked that for these massless embeddings the chiral condensate is non-zero: 

To summarize, the large /i linear analysis has shown that a continuous family of 
massless configurations, parametrized by T < sin~^ {1/2 fi) and carrying a finite R- 
charge density, emerges at the half- integer values of fiR indicated above. In contrast, 
the non-linear DBI equations show a similar set of (unstable) solutions spread out 
in a small window in the vicinity of these critical values. 

The analytical expressions for the mass and the chiral condensate in Eg. ( |5.8| ), 
reproduce quite closely the numerical results displayed in the previous sections (see 
Fig.(0)). 

5.2 Stability analysis: two-point function for ipqipq 

The detailed numerical study of the free energies of the various D7-brane embeddings, 
summarized in Fig.(|D suggests that the values 2/ii? = 3, 5, 7, . . . are critical points 
where two saddle points of the free energy merge. We expect that small fluctuations 
around such points should cost no free energy at linear order, which is essentially what 
our linearized analysis above reveals. If the chemical potential were to be increased 
beyond these critical values, the system should become thermodynamically unstable. 
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Figure 9: Plots of the expressions, obtained in the linearized approximation. 

The natural object to compute, to look for the source of instabilities in the system 
is the two-point function of the operator dual to the slipping mode. In particular, 
the location of the poles in the frequency-space correlator of ipqipq will reveal the 
presence of instabilities in the system. Actually, the linearized analysis above can 
naturally be adapted to allow for a holographic evaluation of the two point function 
of the operator dual to 0. 

Since the dual field theory is formulated on the three-sphere, it is natural to 
consider temporal correlations of the s-wave mode oi O = ipqipq'- 

We begin by expanding the DBI action to quadratic order in 6, allowing for the 
slipping angle to be time dependent, so that 9 = 9{z,t). The induced metric on the 
D7-brane is then 



D7 



+ 299' dt dz + ^— rffi^ + cos^ 9da 



Z'^ 



At quadratic order, the DBI action yields 
^DBI ^ (5.12) 



D7 



dtdz ^ 4^JU + 4^W ^2g/2_ ^ (g2^4^2g2)_3g2^ 



z^ \ (1 + i;z2^2 



This is precisely the quadratic action for a charged scalar in global AdS^ with 
"^scalar = ^3, coupled to a constant, background gauge potential 2/x. This is natural, 
given the interpretation of clS cl chemical potential for a global f/(l) symmetry 
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in the boundary theory. The extra factor of 2 is the charge of the fermion bihnear 
under the U{1) R-symmetry in the boundary. This type of system is also encoun- 
tered in holographic superconductor models, where there is typically a dynamical 
bulk Maxwell field coupled to a charged scalar [11, 23 1, leading to a VEV for the 
latter. The non-trivial embeddings with an R-charged condensate appear to be of a 
similar nature albeit in global AdS spacetime without a black hole. 
Upon Fourier transforming to frequency space, 

/oo 
dte-''^'e{z,t), (5.13) 
'OO 

we note that at this order the effect of a non-zero temporal frequency is a simple 
shift, 

/i^^V + w^ (5.14) 

and 9{z, co) satisfies the linearized equation of motion Eq.(p.l[) after the above replace- 
ment. Therefore, for the holographic calculation of the frequency space propagator 
we can simply import the relevant results from the previous section and make the 
replacement jj, — > The near boundary expansion of 6{z,uj), 

e{z, uj)\^^^ - 0o(a;) z + 62(0;) + |(1 - 4/i2)eo(u;) zMn z + . . . . (5.15) 

can be determined precisely from our linearized analysis above. In particular, we 
already know that regularity at the origin of AdS space leads to a relation between 
62 and Go, from Eq.(|5:8|): 



&2{uj) = &o{uj) X (5.16) 
X 3-(4;.^ + u;^-l){^(i^^^)+V^(i±^^)+ln4 + 27^-l}" . 

Substituting this into the linearized action above, we evaluate the action induced on 
the boundary, 

S^arj = 27r'^^j [ ^^e{e,u)e'{e,-u)\^^^ + ... . (5.17) 

where we have not explicitly written out the counterterms. Differentiating twice with 
respect to Oo(w) which acts as the source for i/jqil^g^cu), results in the frequency space 
correlator. We write the result after re-introducing the radius R of the boundary 
three-sphere and ignoring "contact" terms which are regular in u: 

g{uj) = A/'d7(4^' + uj^- R-^) (1 - |i?Vv^j + (1 + 

(5.: 

Despite the appearance of square roots in this expression, the function only has 
isolated simple poles at 



±v/i?-2(2£+l)2-4/i2, £=1,2,.... (5.19) 
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Figure 10: Simple poles of G{u}) on the real axis, migrate to the imaginary axis as /i is 
increased past 

For /i > some of the poles move to the upper half plane, along the imaginary 

axis, reflecting an instability. 

The /i = functional form of the correlator coincides with the expression deduced 
from the weakly coupled field theory^ in Eq. (|3.8| ). However, the normalizations differ 
between weak and strong couplings. Therefore the correlators do get renormalized 
even in the absence of a chemical potential. 

In the strongly interacting system, the chemical potential has the effect of shifting 
the locations of the poles of the Green's function. With /i 7^ 0, there is a qualitative 
change in the physics in going to strong coupling. The discrete excitations at strong 
coupling consist of mesons/fermion bilinears made from the spherical harmonic states 
of elementary fermions. These composite scalars carry twice the R-charge of the 
fermions, and the poles of the correlator show that they possess twice the energy. 
Importantly, the R-charge chemical potential affects the masses of these states in 
the same way that it would alter the mass of a free complex scalar. A large enough 
chemical potential (bigger than the mass to charge ratio for these states) leads to a 
negative mass squared and triggers an instability favouring a condensate. 

We expect that the instability persists when a fiat space limit is taken, by for- 
mally letting R become arbitrarily large, so that 

g{uj) Mm (w^ + V) In [uj^ + V) , (5.20) 

which has non-analyticities in the upper half plane. 

It is interesting to contrast the above picture with the physics of the free theory 
where the positions of the poles are unaltered by a non- vanishing chemical potential. 
The primary effect there is the Pauli blocking of low lying states, which in turn 
results in an imaginary part for the retarded Green's function in fiat space beyond a 
kinematic threshold l^l > 2yU (Eq.( p.9|) ). 

6. D5 and D3' flavour branes in global AdS^ 

Most aspects of the analysis we performed for probe D7-branes above, carry over 

^In recent work the meson spectrum of probe branes in global AdS was obtained. The poles 
in the correlator we find here correspond to I — I — Q scalar states in that work. 
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to probe D5 and D3-branes whose excitations in AdS^ describe 2+1 and 1+1- 
dimensional defect CFTs coupled to M = 4 SYM at strong coupling ^ ^ 



In the global AdS context, the defect CFT's are formulated on a spatial 5* C S'^ 
and C S^. The eigenvalues of the Dirac operator on S'^ take integer values 
so that 

6eiS^)=iR~\ ^=1,2,... (6.1) 
For the 1+1 dimensional case, fermion energies on the C are half-integer valued 

ee{S') = {i-l)R-\ i = l,2,... (6.2) 

Interestingly, the probe D-branes exhibit the same kind of instabilities we saw above 
for the D7-branes, but now at values of the chemical potential set by the fermion 
energy levels on 5*^ and 5*^. We will not repeat the detailed numerical study of these 
examples, since their basic features turn out to be quite similar to the probe D7- 
brane. Instead, we will briefly explain the analytical results for the lower dimensional 
examples and point out features which are distinct from the four dimensional case. 

6.1 D5-brane probes 

The 2+1 dimensional CFT localized on the intersection of a large-A^ number of 
D3-branes with Nf D5-branes has eight supercharges and an SU{2)h x SU{2)v R- 
symmetry. The defect theory has fundamental hypermultiplets, each containing a 
fermion doublet transforming under SU{2)y. We wish to introduce a chemical poten- 
tial for the U{1) subgroup of this SU{2)v R-symmetry under which the fundamental 



representation fermions are charged. As shown in ||2T| the scalars in the fundamental 
hypermultimplets are uncharged under this. In the near horizon limit of this system 
at large- A^, probe D5-branes in (the Poincare' patch of) AdS^ x can therefore be 
chosen to spin along the internal directions so that R-charge is concentrated in the 
fermionic sector of the fundamental matter fields. 

The same setup, when carried over to global AdS spacetime, involves flavour D5- 
branes with world- volume global ^^54 x S"^. The S'^ C wrapped by the branes 
can be parametrized in terms of the angular coordinate 6 as, 

ds% = de^ + sin^ e dnl + cos^ e d^l . (6.3) 

9 plays the role of the "slipping" angle and the D5-branes can be chosen to spin 
along the azimuthal angle inside the transverse two-sphere S"^. The induced metric 
for the spinning D5-brane with world-volume AdSi^ x 5*^ is 



ds' 



(6.4) 

D5 

^ ^ ' Vsin^^ j dt^ + {zH\zf + 1) ^ + f dnl + cos'ednj. 
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6.1.1 Lineeirized results: 

The D5-brane probe action, in terms of embedding coordinates {Ci}, is 

Sd5 = Nf Td5 J V-det*g , ^fm = Nf Td3 An ^ ^Nj TV, (6.5) 

and the resulting non-linear equations of motion can be numerically solved. Similarly 
to the D7-brane case, we look for embeddings of the slipping mode 9{z) in the 
linearized approximation. This mode maps to a dimension two R-charged operator 
in the boundary 3D CFT, namely the quark (fermion) bilinear ipqipq. The solutions of 
interest are those that get to the origin of AdS. At quadratic order the spinning D5- 
brane action describes a charged scalar in AdS4, with m'^caiar ~ in the presence 
of a constant background gauge potential 2/i,. The linear equation of motion, 

16 + 4^2 + 2^ 16 + z^ + 8z^(l + V) .... 
^ - ' ,(16 -.4) ^ + 2 Oiz) = 0, (6.6) 

subject to the boundary conditions, 

e'{z)\^^, = o, e(z)\^^^ = T, (6.7) 

is solved by 

z sin tan-i f - ^) 

The asymptotic expansion near the boundary then yields the mass and the fermion 
bilinear condensate, 

e{z) ~ ^(0) z + e2z^ + ... (6.9) 

^(0) = T '^^y , ^(2) = -T cos (TT/ii?) (6.10) 

Therefore the linearized embeddings are massless only when, 

fxR = iR-^ £=1,2,3,... (6.11) 

As before, these are critical values of the chemical potential beyond which corre- 
sponding fermion bilinears will become unstable, as we show below. For these finite 
density embeddings the chiral condensate is non-zero, 

(^V'g) = -Mds (-1)' £ e Z. (6.12) 
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6.1.2 Stability of = embedding: 

The calculation of the s-wave 2-point function in frequency space, for the fermion 
bilinear O = ipqipq proceeds exactly as in the 3+1 dimensional theory. The result 
is obtained by taking the ratio of the normalizable and non-normalizable modes in 
( |6.10|) , followed by the replacement /i — )■ 

i W + a;2_ We find, 



g(co) = ^4/^2 + COt(f W + (6-13) 

which has simple poles at 

00 = ±v/(2£)2i?-2 - V, £=1,2,.... (6.14) 

Therefore the onset of instabilities occurs when /ii? > 1, when the poles move to the 
upper half plane. Once again the chemical potential acts as a negative mass squared 
for the scalar excitations at strong coupling 

6.2 D3' probe branes 



The theory on the intersection of the D3-D3' system was studied in detail in . In 
fiat space, the D-brane setup involves two stacks with and Nf D3-branes, spanning 
the and X ^ X ^ X ^ X ^ directions respectively. Each stack yields an A/" = 4 

vector multiplet in four dimensions, coupled to SU{N) x SU{Nf) bifundamental 
degrees of freedom localized at the intersection of the two stacks. 

In the 't Hooft large-A^ limit with Nf fixed, the SU{Nf) vector multiplet de- 
couples since its effective coupling A/ = XNf/N is vanishing in the large- limit. 
Therefore, one obtains Af = 4 SUSY Yang-Mills theory in four dimensions coupled 
to Nf two dimensional A/" = (4, 4) fundamental hypermultiplets localized at the in- 
tersection of the two sets of D3- branes. In this limit, and at strong 't Hooft coupling, 
the embedding of the Nf probe D3'-branes in the AdS^ x geometry induces an 
AdS^ X 5*^ metric on the probes. The C S^, wrapped by the D3' s can be specified 
as usual by the slipping angle 9 

dnl = dO^ + sin^ 6 d^l + cos^ 6 d^l . (6. 15) 

The global isometry of the coordinates transverse to all the D3-branes is 5*0(4) ~ 
SU{2)l X SU{2)ji. Under this isometry, the fermions in the fundamental (4,4) hy- 
permultiplets transform in the (|,0) and (0, |) representations respectively. Their 
scalar superpartners are uncharged under this. The fermion bilinear, which we will 
continue to denote schematically as ipqipq, dual to the slipping mode, transforms in 
the {\,\) representation of SU{2)l x SU{2)r. It has scaling dimension A = 1, 
consistent with a fermion bilinear in two dimensions, and is a |-BPS operator. 

The above system in global AdS spacetime is dual to the 1+1 dimensional defect 
CFT formulated on an equatorial circle of the boundary three-sphere. Imparting an 
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R-charge to the flavour fermions by rotating the D3"s along an angular coordinate 
4> in the with angular velocity the induced metric on the D3' becomes 

ds^ = (6.16) 

D3' 

_ / [L+,Z ) _ ()\dt2^ ^^2^/(^)2 + 1) ^ + f ^ dnj + cos' Od^l 

\ ^ / z z 

This yields the probe D-brane action 

/N 
d^eV-det*^, Td3' = ^. (6.17) 

The slipping mode, at linearized order, is a scalar with ^g^aiar ~ ~^ AdS^. 
It therefore saturates the relevant Breitenlohner-Freedman bound with boundary 
asymptotics, 

e{z) ^(0) zhiz + 6(2) z + ... (6.18) 

Unlike the higher dimensional cases, the term linear in z is now normalizable, whilst 
the non-normalizable mode is the coefficient of the asymptotic logarithmic term. The 
latter plays the role of the hypermultiplet mass. When // » 1, it is consistent to 
linearize the DBI equation of motion, 

(16 + 3z^) 16 + 8^^(1 + 8/.^) + ^^ 

^ - z{lQ - z^) ^ ^2(4 + ^2)2 ^(^) - ■ (6-19) 

The solutions to this equation, which fall in the category of "AdS-filling" cmbed- 
dings, can be written in terms of hypergeometric functions satisfying the boundary 
conditions, 

^(^)L=2 = ^ ^'(^)L=2 = 0- (6-20) 
The corresponding solution is 



(1 + i^ 



1 ^2\2 



1 - i^2)2 • 



4 



From the asymptotic expansion of this function, we read off the fermion mass and 
condensate 

2 

^(0) = -- ^ cos(7r/ii?) , (6.22) 
^(2) = -^T cos{7riiR) [V (i + //i?) + V (i - I^R) + 27£; + ln4] . 
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As for probe D5 and D7-branes, the VEV and the mass satisfy a relation 



0[2) = \ 0(0) (1 + fiR) +iJ - fiR) +2JE + In 4] . (6.23) 

Therefore the embeddings yield zero mass hypermultiplets when 

2/ii? = 2£-l, £ = 1,2,3,..., (6.24) 

which is in agreement with what we would expect for fermion bilinears on 5*^. To 
correctly calculate both one- and two-point functions of the operator ijjqipq, it is useful 
to recall the holographic renormalization procedure for probe branes applied 
to the D3-D3' system. This is particularly relevant in this case due to putative 
logarithmic divergences in 1+1 dimensions. The regularized DBI action is obtained 
by performing a subtraction at z = e ^ 1, 

= 4n'NfTr,-s' J dt (^j\z ^/^d^g - + ^J{e,tf (l + j^)) (6.25) 

Using this, it follows that 

{ifjg^q) = hm = -—^6(2), 6.26 

e-^o e y/—'ydO[e,t) 2ttR 

where we have restored explicit dependence on the radius of the boundary circle. 
Therefore for the massless, finite density embeddings in the vicinity of 2fiR = 2£ — 1 
in the linearized approximation, the chiral condensate is 

6.2.1 Two-point function around = 0: 

It is now straightforward to obtain the two-point function of the fluctuations of 
O = 'ipqipq, about the trivial zero-density embedding in frequency space. A time 
dependent fluctuation 9{z,t) about the 9 = solution, can be expanded near z = 
as (in frequency space), 

e{z,u) = Qo{uj)zlnz + Q2{uj)z + ... (6.28) 

The VEV of 0{u) and the source Qo{^) obey the relation ( |6.23| ) with the replacement 
/i — )■ ^^/Afi"^ + u"^. The Green's function for O is then given by 

(OMO(-.))=-i|^ (6,29) 
2 dK3o[—u) 



A'kR 



^ (^i + |i?VV^^) + ^ (i - \RVW^) + 27i? + ln4 
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As was noted in |^2[) the Green's functions and indeed, the DBI action associated 
to probe D3-branes (|6.17|) is independent of the 'tHooft couphng A, suggesting a 



non-renormahzation theorem. However, this could only possibly be true at /i = 0, in 
the supersymmetric theory. In the situation with a chemical potential we see that 
the locations of poles in the correlator depend non-trivially on fi. In particular, as 
in the higher dimensional examples, the poles of the Green's function shift to 



u = ±v/(2£ - - V, £ = 1,2,3,... (6.30) 

which imply instabilities when /ii? > i. 



7. Discussion 

We have seen that spinning Dp-branes in global AdS^ x S^, dual to strongly cou- 
pled CFT's on spheres with a chemical potential for fermion flavours alone, exhibit 
instabilities (preceded by metastabilities). It appears that in all cases the zero den- 
sity state becomes metastable or unstable when the chemical potential approaches 
integer or half-integer values, depending on the number of spatial dimensions in the 
boundary CFT. At each of these points, the systems would prefer to be in a configu- 
ration with finite R-charge density and with a homogeneous condensate made from 
harmonics of the elementary fermions on the spheres. 

However, the probe DBI analysis leaves the central question unanswered, namely, 
what is the end-point of the unstable directions. The answer to this is not clear from 
the above computations, since the fully non-linear DBI solutions show a runaway 
effective potential. One plausible explanation could be that determining the true 
ground state might require new active "light" modes on the probe brane. This could 
happen, for example, if the ground state were inhomogeneous, as is the case in certain 



lower dimensional models 39 



However, a possible answer may lie in the differing, and somewhat comple- 
mentary, physical pictures that emerge at arbitrarily weak or zero coupling and 
strong coupling. In the free theory, the relevant degrees of freedom are the elemen- 
tary fermions (and their spherical harmonics), whereas in the strongly interacting 
AdS/CFT dual, the D-brane "slipping mode" is the natural degree of freedom. This 
is a composite scalar - a fermion bilinear. A large enough R-symmetry chemical 
potential, effectively induces a negative mass squared for this composite scalar. This 
situation is reminiscent of the Gross-Neveu model where the elementary fermions 
are traded for a composite scalar and an effective potential for this composite scalar 
is subsequently induced, leading to chiral symmetry breaking. In that case, the full 
effective potential requires resumming fermion (flavour) loops. It is conceivable that 
in the present context similar, subleading corrections involving fermion loops need 
to be incorporated into the picture. This would also potentially mean inclusion of 
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back-reaction effects of ffavour branes on tlie background geometry. Tlie possibil- 
ity remains, of approacfiing tfiis question from tlie weakly coupled boundary theory 
where effective four-fermi interactions can be induced via the Yukawa couplings of 
the flavour fermions to the adjoint modes of the J\f — 4 theory on the sphere. 
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Appendix A: Weak coupling Gross-Witten transitions 



In |]T2[, it was shown that large- Yang-Mills theory on 5* with Nf fermions in the 
fundamental representation {Nf/N fixed), undergoes an infinite sequence of Gross- 
Witten transitions as a function of a baryon number chemical potential, at a fixed low 
temperature. In that case, the Euclidean path integral was shown to be dominated by 
a complexified saddle point configuration of eigenvalues of the Polyakov-loop at large- 
N. The complex nature of the saddle point was due to the fact that the Euclidean 
action with a baryon number chemical potential, is not Hermitian - the source of the 
"sign problem". 

Below, we will repeat the essential aspects of this computation for the case where 
a chemical potential is introduced for an axial U{1) symmetry for the fermions. We 
take Nf fiavours of massless Dirac fermions 

{ 1 = 1,2,. ..Nf. (A.l) 

In Lorentzian signature, a chemical potential fi for the axial U{1) symmetry, is 
introduced by 

■t=i 

where the ". . ." represent terms due to any additional matter fields in the theory 
charged under the U{1). To make sense of this we need to assume that the sym- 
metry is anomaly free, which can be achieved in supersymmetric theories with an 
appropriately chosen matter content (e.g. the matter content of A/" = 2 supersym- 
metric gauge theory with SU{N) gauge group and Nf = 2N hypermultiplets) . The 
main feature we want to illustrate is largely independent of these details. 

At finite temperature T, going to Euclidean signature, the theory is formulated 
on S'^ X 5*^. The circumference of the is P = 1/T . Allowing for a Wilson loop 
around the thermal circle, we can choose it to be diagonal and homogeneous on 



1 



Vol(53 



N 

Ao = diag(^i,^2,...M, y"^. = mod27r. (A.3) 



Fermions do not have any zero modes on and at finite temperature. Assuming 
that any other matter fields in the theory do not have zero modes, we can integrate 
out all the Kaluza-Klein harmonics on S"^ to obtain an effective unitary matrix model 
for the Polyakov loop U = expi Aq. In the eigenvalue basis, ( |A.3| ) the effective 
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action for the 9i takes the form pO, p8|, pBl 112 



AT ^ 

E 1^ (1 + e-^^^^-'^)) - (/X ^ -/x) + . . . 



The first term is the Van der Monde pairwise repulsive potential between the eigen- 
values, while the ". . ." represent the contributions from Kaluza-Klein harmonics of 
all other degrees of freedom in the theory. The terms explicitly shown arise from 
integrating out the fermion spherical harmonics, transforming in the fundamental 
representation of SU{N). The energies of the fermion modes and their degeneracies 
on are 

e^ = £+i di = 2i{i+l), £ = 0,1,2... (A.5) 

We note that the effective action is real, in contrast to the situation with baryonic 
chemical potential nevertheless, the physics is quite similar to that case. 

The physics of interest occurs at low T, when /x approaches e^. When /i ~ e^, it 
suffices to focus attention solely on the mode with that energy. In this regime it is 
useful to introduce the fugacity 

Q = e~^''-^^/^ 0«1 for fx<ee, (A.6) 

> 1 for yU > , 

In the large- limit with Nj — )■ oo and Nf/N fixed, the saddle point equation for 
the effective action with /i ~ eg, is 

i=i * 
When fi < €£ , Q ^ 1 and the saddle point equation becomes 

de'p{e') cot (^) ^d,^Q sin^, < 1- (A.8) 

while for fi > Si 

de'p{e') cot {'-^) c^d.-^-sine, o>i. (a.9) 

Here we have introduced the eigenvalue density p{6) on the circle, normalized so that 
J^T, p(^) = 1- These are precisely the equations leading to the Gross- Witten third 
order transition, [RT . 
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This analysis, (from the results of [^) shows that for both fi < ee, {Q ^ 1) and 
/i > Eg, {(i 3> 1), the distribution function p{9) should be "ungapped". 

As /i is increased from low values towards eg, we expect a "gapping" transition 
of third-order when 

deCi^^- (A. 10) 



N 

which corresponds to 



/i = /il = £,-Tln(^^rf,J , (A.ll) 

for small enough T. Therefore for fi > the theory enters a gapped phase for the 
eigenvalues of the Polyakov loop. 

As /i is increased further past Ei, the effective potential becomes of the form 
( [A.9|) and now the theory undergoes yet another transition from the gapped to an 
ungapped phase at 

^(i^-~l A.12 

=^ = fiX = ei + T\nl^di 

Therefore, we learn that the free theory on at very low temperatures, experiences 
a pair of Gross- Witten transitions as a function of the chemical potential fi, when 
the chemical potential is in the vicinity of an energy level of a fermionic mode on 
S\ 

We also expect intuitively that across these transitions, the level with energy eg 
should become occupied. This can be verified by computing the derivative of the 
action/free energy with respect to /x: 

When p < El, Ci <^ 1 at low T, and therefore the corresponding energy level is 
unoccupied. For /i > e^, on the other hand, since 3> 1, the corresponding energy 
level has occupation number, 

N,Ni,Q 1 _|±||^ « N,Nd, (A.14) 

It is also interesting to note that in the ungapped phase, when p{9) ^ the 
Polyakov loop is 

1 

— (Trf/) = / dOe'^pie) ^ , when > 1 or < 1 . (A.15) 

In the gapped phase, however, the Polyakov loop will have an expectation value of 
order one. Therefore, in addition to a jump in the (chiral) fermion occupation number 
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when /il ^ /i ^ the Polyakov loop exhibits a spike across these transitons. The 
situation is illustrated in Fig.(|D. 

Appendix B: Finite density propagators 

In this section we use real time fermion propagators at finite density to sketch the 
perturbative computation of the correlator of the operator il^qipq. We begin by quoting 
the results for fermions in flat space and modify them slightly to account for 
the fact that in our problem the chemical potential is for an axial U{1) symmetry. 
We use the flat space results for the two-point function of ipqipq to infer its form on 
S"^. We package the fermion flavours in the D3-D7 setup into Nf Dirac fermions 
as in Eq.( [A.l| ) {tp^}, i = 1,2, . . . Nf. The momentum space Feynman propagator is 
deflned via 

tS{k;t,0)Sl = [ d'xe-"^-^{ij\x,t)iJmm. (B.l) 



where l,m = l,2...Nf. The Feynman propagator can be written as the sum of 
retarded and advanced pieces 

S{k; t, 0) = S>{k, t, 0) e(t) + S<{k; t, 0) 0(-t) . (B.2) 

For the case with an axial U{1) chemical potential, we write 

5>«)(^;t,0) = i^S>«)(^;t,0) + i±^S>«Hfc;t,0). (B.3) 
With massless, chirality "plus" fermions, = {k,k), we have 

Sl{k-t,Q) = [f^[l~nF{k,^l)]e-'^' + ^,k^,nF{k,^^)e^^') (B.4) 

Stik;t,Q) = ^ + 70^0 [l-nF(A;,/i)]e*'=*) . 

The right-handed fermion propagators 5*^ and are the same as above with the 
replacement /i — )■ — /i. The Fermi-Dirac distributions which turn into step functions 
at zero temperature are given as 

= ^^(fc^),)^^ ^ Q{fi-k) as T^O, (B.5) 
nF{k,fi) = ^^(fc^),)^^ ^ 9(-/i-fc) as T^O. 

We are interested in computing the spatially homogeneous (zero external momentum) 
part of the two-point function for O, 

Nf 

O^^V^g.V^;, Qit)^ I d'x{0{x,t)Om. (B.6) 

i=i 
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Time-ordered propagator: 



It is straightforward to show that 



(2vr)= 



Tr[5(A;;0,t) S{k]t,Q)\ 



(B.7) 



where the trace is over Dirac and colour indices, and all propagators are time ordered. 
On S^, the integral over spatial momenta should be replaced by a discrete sum over 
the spherical harmonic index I with degeneracy 2i{i+l). The Feynman propagator 
for O is then, 



Git) = Nf 



(27r) 



3 T:[5>(A;; 0, t) S<{k; t, 0) e{-t) + S<{k; 0, t) S>{k; t, 0) e{t)] . 

(B.8) 

After some algebra at T = 0, we find, assuming that // > 0, and the number of 
colours N 



g{t) ^ NfN 



d^k 
(2^ 



( [(1 - e(/x - fc)) e^""' + e{ij-k) e-^^'^* ] e(-t) (b.q) 



+ [e(/x - k) e^'^' + (1 - Q{^x - k)) e-^'^' ] Q{t)) . 
The frequency space propagator. 



g{u) 



dte'^'git) 



is therefore, 

+e{li - k) 



d^k 
1 



u -\-2k — ie uj — 2k + ie 

1 1 

+ 



+ 



(B.IO) 



(B.ll) 



^uj — 2k — ie uj + 2k — ie u — 2k + ie uj + 2k + ie^ 
Rewriting the energy denominators using l/{x — ie) = P(l/x) + i7:S{x), we get 



g{uj) ^ NfN 



_fk_ 
(2^ 



uj + 2k uj-2k 



(1 - 2e(// - k)) 



(B.12) 



Recalling that k = \k\, the spatial momentum integrals can be performed so that. 



g{u) = 



167r2 



a; Incj +0; In 



cj^ - 4//' 



2\ 2 



(B.13) 
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Here we have ignored additive constants and contact terms, analytic in oj. The chem- 
ical potential has modified the /i = result for the Green's function of introducing 
branch points at w = ±2yU. 



Retarded propagator: 

From the expression for the time ordered propagator, we see that the contribu- 
tions of on-shell intermediate states, even in the presence of the chemical potential, 
only yield singularities on the real u axis. Therefore the different propagators - 
time-ordered and retarded, differ only in the choice of ie prescriptions for these sin- 
gularities. We confirm this by carrying out the explict computation of the retarded 
propagator, defined as 



We find that 
Gr{uj) = NfN 



ij dte'^^'Qit) / c/^x([C>(f,t),C(0)]). 



(B.14) 



(2^ 



1 



1 



uj + 2k + it tu — 2k + ie 



{1 - 2e{fi - k)) . (B.15) 



The real part of this expression (for real u) is exactly the same as the time-ordered 
propagator. The main difference lies in the imaginary part which is determined by 
the ie prescriptions. Up to additive constants, the retarded propagator at finite 
density in fiat space is 



NfN 
' 167r2 



w In w + In 



2\ 2 



+ 2i7iujhgn{uj)e{uj^ - 4/i2)(B.16) 



We can also deduce the form of the two-point function of O on from the 
formula ( [B.12D by replacing the integrals by a discrete sum over the energy levels on 



NfN 



1=1 \ < t- 



1 



tu — 2ef 



2eifi - Si) 



u + 2ee tu — 2et 



where for fermion harmonics 

ei = {l + \)R-\ de = 2i{i+l) 
Appendix C: D7-brane embedding ansatz 



(B.17) 
(B.18) 



In this appendix, we see why the embeddings with g{z) ^ are not suitable for the 
theory with massless hypermultiplets. With the general ansatz 



(l){z,t) = 2fit + g{z) 



(C.l) 
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the DBI action for the probe D7-brane becomes 



(I _ 

Sm = 2ix^Md7 I dt dz cos^ 6 ^ x (C.2) 

z^ 



1 + zW{zy) (^(1 + - V ^2 sin^ 9^ + g'{zy sin^ 9 z^l + 
As usual, since the action does not depend on g{z)., we have 



^ - = 27T MdtC = constant. (C.3) 
dg' 



From this condition we find, 



1 + zW{zy) (^{1 + 4)2 - V ^2 sin2 ^^ 
\| cos6esin2^(l + 4)'(1 - t)^^"^ sin^(l + 4)^ 



(C.4) 



The equation of motion for 9{z) is best obtained through the Legendre trans- 
formed action 

S..^S..-2.^M..I,U.,'U)c. (C,5) 

which simphfies to 



S'lt = 271^07 dtdz X (C.6) 

J 2sin^(l + ^) 

'l + zW{zy) (^(1 + 4)^- V^^ sin^^) (cos6 ^sin2 ^(1 + 4)^(1- 4)® ^"^-^2) 

From the equation of motion for 9{z), we find that the asymptotic behaviour of the 
shpping angle near the boundary is 

9{z) = ^(0) z + 9i2) + i^(o) (1 - V) ^Mn 2 + . . . , (C.7) 

so that the constant of the motion c does not appear in the mass or the condensate 
of the corresponding field theory operator. However, it does make an appearance in 
g{z) and therefore in (f){z,t), which then has the asymptotic form. 



(f){z,t)^2fit + --^z + ... (C.8) 

'^(0) 

For massless fields this implies a putative divergence in (p, near the boundary. This 
suggests that non-zero values of c are potentially problematic for configurations 



with massless hypermultiplets and so c must be set to zero for this case. (See 
for a detailed explanation of related issues.) 
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